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Abstract. The first phase of the Alborz Observatory Array (Alborz-I) consists of 20 plastic 
scintillation detectors each one with surface area of 0.25 m 2 spread over an area of 40x40 m 2 
realized to the study of Extensive Air Showers around the knee at the Sharif University of 
Technology campus. The first stage of the project including construction and operation of a 
prototype system has now been completed and the electronics that will be used in the array 
instrument has been tested under held conditions. In order to achieve a realistic estimate 
of the array performance, a large number of simulated CORSIKA [1] showers have been 
used. In the present work, theoretical results obtained in the study of different array layouts 
and trigger conditions are described. Using Monte Carlo simulations of showers the rate 
of detected events per day and the trigger probability functions, i.e., the probability for an 
extensive air shower to trigger a ground based array as a function of the shower core distance 
to the center of array are presented for energies above 1 TeV and zenith angles up to 60°. 
Moreover, the angular resolution of the Alborz-I array is obtained. 
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1 Introduction 

Study of the primary cosmic rays provides important information on their origins and accel¬ 
eration mechanisms. The energy spectrum of cosmic rays follows a power law and shows a 
bend around 3xl0 15 eV, known as the knee region, where the spectral index changes from 
-2.7 to approximately -3.1 [2], At these energies, indirect measurements of cosmic rays due 
to their low flux are performed using ground based arrays of particle detectors. In the last 
decade, findings of many high-profile experiments confirm the general view which attributes 
the knee to processes of magnetic confinement occurring either at acceleration regions or as 
diffusive leakage from the Galaxy or both. 

Air shower measurements can determine the arrival direction, energy, and mass of the 
primary cosmic ray. A particle detector array provides statistical samples of secondary parti¬ 
cles in the shower disc, and the resolution achieved in the shower reconstruction depends on 
the sampling fraction of secondary particles from the shower disc [3]. 

Lowest and highest energies at which cosmic rays can be detected by arrays depend on 
the distance between neighboring detectors and overall size of detector array respectively. 
Therefore geometry and design of an array have an important role to extend the detectable 
energy range. In this paper, geometrical effects of array layout on the number of detected 
events, trigger probability and angular resolution are studied by increasing distance between 
neighboring detectors and also changing configuration of the array. 

In order to show dependence of the trigger probability on the primary particle energy 
(E) and shower core distance from the array center (r), detailed studies on the trigger rate for 
two rectangular and cluster layouts are carried out. The trigger probability of the simulated 
array as a function of r and E, known as the trigger probability function P(r, E), is presented 
as a fundamental parameter to compare different trigger conditions in the both rectangular 
and cluster layouts. On the other hand, the angular resolution of the array for different 
trigger conditions and layouts is calculated using a chi-squared minimization algorithm. Two 
different shower sets simulated by CORSIIvA, which differ in the distribution of energy and 
zenith angles, have been applied for these studies. 
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Figure 1 . Cluster layout of the Alborz-I array where 20 SDs will be deployed in a square area. 
Distance between neighboring detectors in the central pentagon and also in the surrounding triangles 
around the central pentagon is 5 m. 


2 The Alborz-I experiment 

The Alborz Observatory Array (AOA) has been designed to measure the flux, arrival direc¬ 
tion and energy of the primary cosmic ray around the knee. The AOA has been planned to 
have both Scintillation Detectors (SDs) and Water Cherenkov Detectors (WCDs). 

The Alborz-I as the first phase of the AOA is realized at Sharif University of Technology 
(SUT), Tehran (35°43 , N 51°20 / E, 1200 m a.s.l corresponding to an average atmospheric depth 
of 890 gem -2 ), and consists of 20 scintillation detectors (each one 0.25 nr 2 ) covering an area 
of about 1600 m 2 . In the next phase, 10 WCDs will be added to it. The used detectors are 
optimized by a series of experiments and simulations [4, 5]. 

Each scintillation detector of the Alborz-I has a pyramidal galvanized iron light enclo¬ 
sure with 1 mm wall thickness and diffuse white coated inner surface, housing a 50x50x2 
cm 3 NE102A scintillator. Each detector is embedded in a station with a 5 mm cement ceiling. 
Scintillation light emitted due to passing a particle through the scintillator is read out from 
tip of the pyramidal light enclosure by a 5 cm photomultiplier tube (PMT, 9813B) at the 
distance of 20 cm of the scintillator sheet for timing measurements. 

Two different array configurations, one cluster and other rectangular grid, are proposed. 
As shown in figure 1, in the cluster layout, 5 out of 20 detectors arranged on an internal 
pentagon with side length 5 m and the rest set on 5 surrounding triangles with the same side 
length 5 m which are placed on the corners of an external pentagon with side length about 
18 m. 

In selection chain of the array properties, such as array layout, trigger condition and 
energy range of interest, several parameters including the trigger probability function, the 
number of recorded events and the angular resolution of the array are studied for different 
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implemented conditions in the simulations. The results from study of the array properties 
are mentioned in the following sections. 

3 Monte Carlo simulations 

A realistic estimate of the array response and angular resolution, as two main characteristics 
of an array, requires a detailed knowledge of the physics of shower development and interac¬ 
tion of secondary particles with the detector material. 

Two different EAS sets simulated by CORSIKA (version 6.9), which differ in the distri¬ 
bution of energy and zenith angles, have been used. QGSJET-II [ 6 ] and GHEISHA (E ^80 
GeV) [7] models have been employed for high and low energy hadronic interactions, respec¬ 
tively. At high energies, the showers have been generated in the High Performance Computing 
Center (HPCC) of SUT with respect to time consuming simulations. 

In order to reliable estimate of the trigger probability function, a shower set consists 
of 3600 extensive air showers with a composition of 88 % proton and 12% alpha as primary 
particles have been simulated. In this set, air showers were simulated with azimuth angles 
from 0° to 360° and zenith angles between 0° and 60° distributed as 7oc sin# cos# d# (which 
the sine term respects the solid angle element of the sky, while the cosine term takes into 
account the geometrical efficiency of a flat horizontal detector). The energy of primaries is 
discretely distributed in a range between 10 12 eV and 10 16 eV in steps of 0.5 in log E. For 
the both layouts, the core position of any shower is uniformly distributed at all over the 
array surface. The surface of the array is divided into 25 square pixels, each one with sur¬ 
face area of 49 m 2 , and for any shower the core position is set in the center of each pixel 
in order to increase the statistics of simulated EASs and also to cover different zones of the 
array. Also a wide square grid including maximum number of 2809 square pixels is consid¬ 
ered to study showers arriving inside and outside of the array surface while fired the detectors. 

Another set of showers includes 12000 proton and alpha showers, with the same previ¬ 
ous composition, distributed over a continuous energy range from 2xl0 14 eV to 4xl0 14 eV 
according to a power law with the spectral index of 7 — -2.7 have been generated to estimate 
the angular resolution of the array. This energy range is selected according to the results of 
the first simulated showers set which reveal a significant number of recorded events in this 
range in comparison with other bins of energy. In this set, the arrival direction of showers 
is selected from 0° to 60° for zenith angle in steps of 5° and from 0° to 360° for azimuth angle. 

The response of the detectors to gammas as secondary particles is calculated by means of 
GEANT4 detector simulation toolkit [ 8 ]. The energy deposition of gammas in the scintillation 
detectors is computed using simulated air showers generated by CORSIKA which are used as 
input for detailed GEANT4 simulations. 

4 Trigger probability function 

The trigger probability of a ground based array depends on several independent physical 
parameters: i) the characteristics of the primary cosmic ray which initiates an air shower, 
e.g., energy and mass of the primary, ii) the type of the array detectors, iii) the trigger 
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condition used to detect air showers, iv) the array layout, v) the geometry of the incoming 
shower, e.g., its incidence zenith angle and distance of shower core to the array center [9]. To 
show these dependencies, the trigger probability function, P(r,E), has been presented as a 
fraction of showers in the energies of interest (E) which fulfill the trigger condition in different 
bins of shower core distance from the array center (r), as 

r ^ _ N trigger jr, E) 

N'incident (j" > 

Where Nt r i gger (r, E) and Ni nc id en t{r, E) are number of showers with energy E which 
fulfill the trigger condition at distance r and total number of incident showers with given 
parameters r and E, respectively. 

The trigger probability function is an increasing function of energy, as the shower size 
gets bigger, and a decreasing function of distance from the shower core, as the density of 
particles gets smaller [10]. The trigger probability function changes between 0 and 1, which 
0 happens for large shower core distances and 1 happens for the shower core distances close 
to the array center in high energy showers. 

In order to optimize the array size in a rectangular configuration, four different array 
sizes with detector spacings of 1.5 m, 3.5 m, 7 m and 14 m are studied. It should be noted 
that the surface area of the site at the SUT campus lets the maximum detector spacing of 7 
m. However, in the last phase site of the array will be extended. The trigger conditions are 
provided by randomly at least 4-fold (4 out of 20), 6-fold (6 out of 20), 10-fold (10 out of 20) 
and finally fully triggering detectors. Due to the small surface of detectors, the fully triggering 
condition can provide more satisfactory statistical sample of charged particles in the shower 
disc for future studies on the energy resolutions and mass composition. Although for the 
fully triggering condition, the statistic of recorded events decreases significantly. The trigger 
conditions at least 4-fold and 6-fold achieve larger statistics of recorded events. Though, as 
discussed in the next section for the at least 10-fold condition, the angular resolution accuracy 
improves significantly in all zenith angles rather than at least 4-fold and 6-fold conditions. 
Thus, in the rectangular layout, we focus the analysis on the at least 10-fold condition. 

The trigger probability function from simulations for different given detector spacing in 
different bins of energy for the rectangular layout is shown in figures 2 and 3. As can be 
seen, the maximum shower core distance ( r max ) for detection in each energy is approximately 
independent of the detector spacing while it increases with energy for a given detector spacing. 

Also the minimum energy of a shower that can be detected by the array obviously de¬ 
pends on the detector spacing so that for the smaller distances between neighboring detectors 
in the rectangular grid, efficiency on lower energies increases significantly. In other words, 
the smaller array sizes allow to detect lower energies as the shower disc size becomes smaller. 
Therefore to focus on the knee, array size should be large enough. 

Beside the rectangular grid, a cluster layout is also studied. Through several imple¬ 
mented trigger conditions, those conditions with full triggered central cluster and two sur¬ 
rounding triangles and or three surrounding triangles and or four surrounding triangles, (i.e., 
5 out of 5 in the central cluster and 2 adjacent triangles, 5 out of 5 in the central cluster and 3 
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Figure 2. Array trigger probability as a function of shower core distance to the array center (r) 
for the 10-fold condition at different energies, from 3xl0 13 eV up to 3xl0 15 eV in steps of 0.5 in the 
logarithmic scale, for detector spacing of 1.5 m and 3.5 m in the rectangular grid. 
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Figure 3. Array trigger probability as a function of shower core distance to the array center (r) 
for the 10-fold condition at different energies, from 10 14 eV up to 3xl0 15 eV in steps of 0.5 in the 
logarithmic scale, for detector spacing of 7 m and 14 m in the rectangular grid. 


adjacent triangles, and finally 5 out of 5 in the central cluster and 4 adjacent triangles) show 
higher efficiencies. In addition, their trigger probability are approximately the same. 

Therefore, in triggering selection, these three trigger conditions are selected: randomly 
10-fold, 5 out of 5 in the central cluster, and finally 5 out of 5 in the central cluster and 
3 adjacent triangles. As shown in figure 4, two last trigger conditions in the cluster layout 
decrease the maximum shower core distance ( r max ) at different energies, e.g., r max is limited 
to 100 m for the primary energy of 3xl0 15 eV by triggering 5 out of 5 in the central clus¬ 
ter. Also the condition of randomly at least 10 -fold triggering shows almost the same trigger 
probability function in the cluster layout as well as the rectangular grid with spacing 7 m 
between neighboring detectors. 

As shown in figures 3 and 4, the trigger probability function becomes almost 0.8 at the 
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Figure 4. Array trigger probability as a function of shower core distance to the array center (r) at 
different energies for triggering 10 out of 20 (top), 5 out of 5 in the central cluster (left-down), 5 out 
of 5 in the central cluster and 3 adjacent triangles (right-down). 


energy of 3xl0 15 eV for shower core distances close to the array center for detector spacing 
of 7 m and 14 m in the rectangular grid as well as the cluster layout. 

In addition to the trigger probability function, number of triggered events per day in 
different bins of energy ( N(E , E + A E)) that fulfill the implemented trigger condition for 
different array layouts can be calculated by a convolution of the trigger probability function 
P(r, E) and differential flux of primary cosmic rays <P(E) (per area, time, solid angle and 
energy) (eq.(4.2)), 


N(E, E + AE) 


E+AE 


P(r,E)$(E) dS dE dn dt 


S JE 


(4.2) 


where the intervals are selected according to the default initial conditions in the simu¬ 
lations. 5 is the surface of a square grid which considered to study showers arriving inside 
and outside of the array while triggered it. <P(E ) (eqations (4.3a), (4.3b) and (4.3c)) is a 
primary all-particle spectrum in the knee region determined with the Tibet air shower array 
which is located at an almost ideal atmospheric depth for this energy range and is highly 
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Detector spacing 3.5 m 




Figure 5. Energy distribution of triggered events per day for different spacing between the neigh¬ 
boring detectors in the rectangular layout for the 10-fold condition. 


instrumented [10]. In the energy range below the knee, for E< 5.62-10 14 eV , 

HE) = 1.5 • 10 ~ 2 °( 5 )~ 2 - 60± °'° 4 [rn-h-'sr-'eV- 1 ], (4.3a) 

and above the knee, for E > 7.08 • 10 15 eV, 

HE) = 1.2 • 10 _23 (-3.00±0.05 J 2 -1 -1 y- 1 ] ( 4 . 3b ) 

v > v 7.08-10 15; l \1 \ ) 

and in the knee, at E = 1.78 • 10 15 eV, 

HE) = 6.7 • 10" 22 [m - 2 s - 1 sr - 1 eV -1 ]. (4.3c) 

Figure 5 shows the energy distribution of triggered events per day that fulfill the 10-fold 
condition for different array sizes in the rectangular layout. The results show an energy de¬ 
pendence of maximum triggered events on the detector spacing so that it changes from 10 14 
eV in the small array sizes (detector spacing of 1.5 nr and 3.5 m) to 3xl0 14 eV in the large 
enough array sizes (detector spacing of 7 m and 14 m). 
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Figure 6. Energy distribution of triggered events per day for different trigger conditions in the 
cluster layout. 


The energy distribution of triggered events per day for different trigger conditions in 
the cluster layout is shown in figure 6. As can be seen, the maximum happens in 3xl0 14 
eV so that the array size is the same as the rectangular grid with detector spacing 7 m. It 
should be noted that the errors are due to systematic uncertainties of the cosmic ray spectrum 
(eqations (4.3a), (4.3b) and (4.3c)) which is used in the Monte Carlo simulations to estimate 
the number of triggered events. 


5 Angular resolution 

The angular resolution of the Alborz-I is discussed by considering the second shower set in¬ 
cluding 12000 proton and alpha showers (with the same previous composition) distributed 
over a continuous range from 2xl0 14 eV to 4xl0 14 eV with 6 between 0° and 60° in steps 
of 5°. This energy range is selected according to receiving maximum number of events in 
E 1 —3xl0 14 eV. 

Using a chi-squared minimization algorithm the angular resolution of the array for dif¬ 
ferent trigger conditions and layouts is calculated. The angular resolution is directly related 
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Figure 7. Azimuth angle resolution for the detector spacing 3.5 m (top), 7 m (left-down) and 14 m 
(right-down) in the rectangular layout. 


to the accuracy of timing measurement. Therefore to provide more realistic results in the 
simulation, the uncertainty on the timing measurement is implemented and also core posi¬ 
tions are randomly distributed at all over the array surface. 

The timing measurement uncertainty by the Alborz-I instruments is calculated by: 

®instr — \/'^det.^ T ^ elec. ^ 1.82 US, (5.1) 

where cr e ; ec . = 0.10 ns is the electronic uncertainty and (Jde.t. = 1-82 ns is the detector 
uncertainty [5]. 

As shown in figures 7 and 8, the angular resolution is obviously improved for triggering 
10 -fold in comparison with two other trigger conditions for different array sizes in the rectan¬ 
gular grid. For the detector spacing 1.5 m due to small size of the array with respect to the 
size of simulated shower disc, the angular resolution is minimized and not shown here. 

Figure 9 shows that the angular resolution is an increasing function of the array size 
and although the angular resolution improves for the detector spacing 14 m, but there is a 
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Figure 8. Zenith angle resolution for detector spacing 3.5 m (top), 7 m (left-down) and 14 m 
(right-down) in the rectangular layout. 


small difference between 7 m and 14 m detector spacing specially for zenith angles greater 
than 10°. In addition as mentioned above, the site of array limits detector spacing to 7 m. 

As illustrated in figure 10, the angular resolution improves in the pentagon layout in com¬ 
parison with the rectangular grid with triggering 10-fold and detector spacing 7 m. Through 
the three selected trigger conditions in the pentagon layout, those which include fully triggered 
central cluster show almost the same angular resolutions. The errors represent statistical un¬ 
certainties as well as systematic uncertainties which is due to the time resolution of the array 
instruments. 


6 Conclusion 

Monte Carlo studies based on the CORSIKA program on the Alborz-I array characteristics 
have been performed. The trigger probability function, angular resolution and energy distri¬ 
bution of recorded events per day have been estimated for two layouts and different trigger 
conditions. 
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Figure 9. Azimuth angle resolution (left) and zenith angle resolution (right) as a function of 
primary’s zenith angle in the rectangular layout for the 10-fold condition. 
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Figure 10. Azimuth angle resolution (left) and zenith angle resolution (right) as a function of 
primary’s zenith angle for the pentagon layout. 

The array response, trigger conditions and shower development are summarized in the 
trigger probability function. The results indicate that although at lower energies the trigger 
probability decreases significantly but due to high flux of cosmic rays, large number of events 
are recorded. Moreover at higher energies (around the knee of the cosmic ray spectrum), the 
array is more efficient. 

In the rectangular layout, four different trigger conditions and four array sizes were 
considered. Study on the angular resolution showed that it obviously improves as the array 
size and number of triggered detectors increase. But limitation in the surface area of the 
array and, the significant decrease in recording events for fully triggering of the array ( 20 -fold 
triggering condition) made us to consider the detector spacing 7 m and the triggering of at 
least 10 -fold as a default condition to compare with the results of the pentagon layout. 

In the pentagon layout the implemented trigger conditions, which include fully triggered 
central cluster together with some of the adjacent triangles yield approximately the same 
efficiency in different energy ranges. Energy distribution of recorded events for different trigger 
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conditions in the pentagon layout maximize at around 3xl0 14 eV. The angular resolution of 
the pentagon layout improves significantly in comparison with the rectangular grid with 7 m 
detector spacing and 10 -fold triggering conditon. 
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